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Abstract 

We search for the decay B° — > p°p° and other possible charmless modes with a 7r + 7r~7r + 7r~ 
final state, including B° -> p°/o(980), B° -» f (980) f (980), B° -» / (980)^vr, B° -> p°Tnr and 
non-resonant B° — > Air. These results are obtained from a data sample containing 520 x 10 6 
BB pairs collected by the Belle detector at the KEKB asymmetric-energy e + e~ collider. We 
measure a branching fraction of (0.9 ± 0.4±g;|) x 10" 6 , or B(B° -> p°p°) < 1.6 x 10~ 6 at the 
90% confidence level. The significance including systematic uncertainties is 1.8<r. These values 
correspond to the final state being longitudinally polarized. We also measure the branching 
fraction of non-resonant B° — > Att decay to be (10.2 ± 4.7lj'|) x 10 -6 with 2. lex significance, 
and set the 90% confidence level upper limit B(B° — > Att) < 17.3 x 10 -6 . For the other related 
decays, B° -> p°/o(980), B° -> f (980) / (980), B° -» / (980)vr7r and B° -> /vrvr, no significant 
signals are observed and upper limits on the branching fractions are set. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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INTRODUCTION 



In the Standard Model (SM), CP violation in the weak interaction can be measured 
through the differences between B and B mesons decays. Measurements of CP violation 
help determine (or constrain) the elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix [lj and thus test the Standard Model. The time-dependent CP 
asymmetry for the decay of a neutral B meson via a b —>■ u process into a CP eigenstate 
can determine the CKM phase angel 02 = ~ v td v* ■ However, in addition to the b — > u 



tree amplitude, a b — > d penguin amplitude contributions, and thus an isospin analysis [2 
is needed to determine 02- 

At present, the 2 constraints on B — > ttk pjj, B — > pn [4] and B —>■ pp [5] are well 
studied. However, it is necessary to measure the branching fraction and polarization 
of B° — > p°p° for improved 02 constraints. Angular analysis can provide additional 
information on VV decays such as B — > pp. Polarization measurements in the B° — > 
p + p~ and B + — > p + p° modes [5| show the dominance of longitudinal polarization, thus 
B° — > p + p~ is a CP eigenstate; measurements of the branching fraction and polarization 
of CP- violating asymmetry in B° — > p°p° decays would complete the isospin triangle and 
improve the constraints on 2 . 

Theoretically, the tree contribution to B° — > p°p° is color-suppressed and thus its 
branching fraction is much smaller than that of B° — > p + p~ or B + — > p + p°. The decay 
rate for p°p° is sensitive to the penguin amplitude. Predictions for B° — > p°p° using 
perturbative QCD (pQCD) [f| or QCD factorization 0, [1] approaches suggest that the 
branching fraction B(B° — > p°p°) is at or below 1 x 10~ 6 and that the longitudinal 
polarization fraction, Jl, is around 0.85. A non-zero branching fraction for B° — > p°p° 
was first reported by the BaBar collaboration they measured a branching fraction of 
B{B° -> p°p°) = (1-07 ±0.33 ±0.19) x 10~ 6 with a significance of 3.5 standard deviations 
(a), and a longitudinal polarization fraction, = 0.87 ± 0.13 ± 0.04. 

A theoretical prediction for the non-resonant B° — > Att branching fraction is around 



1 x 10 4 lOj. The most recent measurement of this decay was made by the DELPHI 
collaboration [ll|, who set a 90% confidence level upper limit on the branching fraction 
of2.3xl0- 4 . 

In this paper, we report the results of a search for B° — > p°p° along with other modes, 
including B° p°/o> B° — > /o/o? B° — > foinc, B° — * p°7nr and non-resonant I? — > An. 

DATA SET AND APPARATUS 

The data sample used contains 520 x 10 6 BB pairs collected with the Belle detector 



at the KEKB asymmetric-energy e + e (3.5 and 8 GeV) collider [12(], operating at the 



T(4S') resonance. The Belle detector [13l. Il4j is a large-solid-angle magnetic spectrometer 



that consists of a silicon vertex detector, a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time- 
of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that provides a 1.5 
T magnetic field. An iron flux-return located outside of the coil is instrumented to detect 
K° L mesons and to identify muons. 

We study signal and backgrounds using Monte Carlo (MC) simulation. For these 
simulations, signal decays, generic b — > c decays and charmless rare B decays are generated 
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with the EVTGEN 15j event generator. Signal MC event gene ration utilizes the PHOTOS 
simulation package to take account of final-state radiation [16]. The continuum MC events 
are generated through e + e~ — > 7* — > qq (q = u, d, s, 
GEANT3 (lH package is used for detector simulation. 



decays in JETSET [17]]. The 



EVENT SELECTION AND RECONSTRUCTION 

B° meson candidates are reconstructed from neutral combinations of four charged 
pions. Charged track candidates are required to have a distance-of-closest-approach to 
the interaction point (IP) of less than 2 cm in the beam direction (z-axis) and less than 0.1 
cm in the transverse plane; they are also required to have a transverse momentum > 0.1 
GeV/c in the laboratory frame. Charged pions are identified using particle identification 
(PID) information obtained from the the CDC (dE/dx), the ACC and the TOF. We 
distinguish charged kaons and pions using a likelihood ratio TZpib = £k/(£k + 
where C^{Ck) is a likelihood value for the pion (kaon) hypothesis. We require TZ-piu < 0.4 
for the four charged pions. The pion identification efficiency is 90%, and 12% of kaons are 
misidentified as pions. Charged particles positively identified as an electron or a muon 
are removed. 

To veto B — > D^n backgrounds, we remove candidates that satisfy any one of the fol- 
lowing conditions: IMO^tt^-M^J < 13MeV/c 2 or IM^tt^-M^I < 13MeV/c 2 , 
where h is either a pion or a kaon, and and Mpo are the nominal masses of D( s ) and 

D°, respectively. Furthermore, to reduce the B° — > a^n T feeddown in the signal region 
(0.55 GeV/c 2 < Mi^^n) < 1.35 GeV/c 2 ), we require that the highest momentum pion 
have a momentum in the T(4S') center-of-mass (CM) frame within the range 1.30-2.65 
GeV/c. 

The signal event candidates are characterized by two kinematic variables: beam-energy 
constrained mass, M bc = a/ -^bcam — ^b> an d energy difference, AE = E B — E hcSLlQ , where 
-^beam is the run-dependent beam energy, P B and E B are the momentum and energy 
of the B candidate in the T(4S) CM frame. We select candidate events in the region 
5.24 GeV/c 2 < M bc < 5.30 GeV/c 2 and |A£| < 0.01 GeV. 

The invariant masses Mi(tt + tt~) and M 2 (7r + 7r~), are used to distinguish different modes 
from signal. In B° — > p°p° — > (7r + 7r - )(7r + 7r~) decays, there are two possible combi- 
nations: (Tri^rX 71 ^ 71 ^) and {^t'^2){' K t' K i)i where the subscripts label the momentum 
ordering (e.g., tt± has higher momentum than )■ According to a MC study, assuming 
B° — > p°p° decays are longitudinally polarized, the reconstructed 7r + 7r _ pair contain- 
ing a high momentum 7r + (7r~) and a low momentum 7r~(7r + ) will corresponds to the 
correct combination 84.8% of the time. Here we consider both possible 7r + 7r~ combi- 
nations and select candidate events if either one of the combinations lies in the p°p° 
signal mass window, which is 0.55 GeV/c 2 < M(n^n^) H M^^n^) < 1-35 GeV/c 2 or 
0.55 GeV/c 2 < M(n^ n 2 ) n M^vrf ) < 1.35 GeV/c 2 . If a candidate event has two tt+tt" 
pair combinations that both lie in the p°p° signal mass window, we cannot distinguish 
which p°p° mass combination is correct. In such cases, we select the 7t + tt~ pair containing 
a high momentum 7r + (7r~) and a low momentum tt~(tt + ) as the correct combination; with 
this selection, 0.2% of the signal is incorrectly reconstructed according to the MC. For 
fitting, we randomly assign the ir + ir~ pairs to either M 1 (7r + 7r _ ) or A-f 2 (7r + 7r~) to sym- 
metrize the 2-D invariant mass distribution. Therefore, the probability density functions 
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(PDF) for M 1 (7r + 7r-)-M 2 (7r + 7r-) are symmetric in both the M x (n + n-) and M 2 (n+n-) 
projections. 



BACKGROUND SUPPRESSION 

The dominant background is the continuum. To distinguish signal from the jet-like 



continuum background, we use modified Fox- Wolfram moments [19j, which are combined 
into a Fisher discriminant. This discriminant is combined with PDFs for the cosine of the 
B flight direction in the CM and the distance in the z-direction between two B mesons 
to form a likelihood ratio 1Z = C S /(C S + C qq ). Here, C s (£ qq ) is a likelihood function for 
signal (continuum) events that is obtained from the signal MC simulation (events in the 
sideband region Mt, c < 5.26 GeV/c 2 j. We also use the flavor tagging quality variable r 
provided by a tagging algorithm [20J] that identifies the flavor of the accompanying B° 
meson m the T(AS) -> B°B . The variable r ranges from r = (no flavor discrimination) 
to r = 1 (unambiguous flavor assignment), and is used to divide the data sample into 
six r bins. Since the discrimination between signal and continuum events depends on 
the r-bin, we impose different requirements on TZ for each r-bin. We determine the TZ 



requirement so that it maximizes the figure-of-merit N s / ^ N s + N qq , where N s (N qq ) is 
the expected number of signal (continuum) events in the signal region (|Ai£| < 0.05 GeV 
and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 ). 

After applying all selection criteria, 17% of selected events have multiple B° — ► p°p° 
candidates. For these events we select a single candidate that having the smallest \ 2 
value of the B° decay vertex reconstruction. The detection efficiency for the signal MC 
is calculated to be 7.11% (9.57%) for longitudinal (transverse) polarization. 



ANALYSIS PROCEDURE 

Since there are large overlaps between B° — > p°p° and other signal decay modes in the 
Mi(nn)-M 2 (nn) distribution, it is better to distinguish these modes using a simultaneous 
fit to a large Mi(nn)-M 2 (nn) region. However, the correlations between (AE, M bc ) and 
(Mi, M 2 ) for backgrounds and signals will lead to large uncertainties as the region in 
M\(nn)-M 2 (nn) used in a simultaneous fit increases. The Mi(nn)-M 2 (nn) distribution 
is shown in Fig. [1] for Monte Carlo samples of non-resonant B° — ► An, B° — > p°nn, 
B° — ► p°p° decays, and the data. The analysis proceeds in three setps: we first measure 
the non-resonant B° — > Atc branching fraction in area A, then measure the B° — ► p°7T7r 
in area B by fixing the branching fraction of non-resonant B° — ► An determined in area 
A. Finally, we determine B° — > p°p° and the other possible decay modes in area C by 
fixing the branching fractions of the non-resonant B° — > An and B° — > p mi decays. This 
procedure minimizes bias caused by unknown correlations between variables (AE, M^ c ) 
and variables (Ml, M 2 ). 

For the non- resonant B° — > Air branching fraction measurement, we use the nomi- 
nal 7T+7T- mass region 1.05 GeV/c 2 < M m (mi) < 1.70 GeV/c 2 and 1.35 GeV/c 2 < 
j ^2(i)( 7r7r ) < 1-70 GeV/c 2 (area A), in which only B° — > a^7r T and non-resonant 
B° — > An signals having the same final state need to be considered. For the B° — > p°nn 
branching fraction measurement, we use events in the mass regions 1.30 GeV/c 2 < 
M m (nn) < 1.70 GeV/c 2 and 0.55 GeV/c 2 < M 2{l) (nn) < 0.95 GeV/c 2 (area B), in 
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which B° — > afn T , non-resonant B° — > Air and 5° — > p°7r7r decays are included. The 
yields for 5° — ► p°p°, p°/oj /o/o an d /o 71 " 71 " decays are measured using a simultaneous fit 
in which the branching fractions of non-resonant B° — > 47T and P° — > p°7T7r in the mass 
region 0.55 GeV/c 2 < Mi i2 (7T7r) < 1.35 GeV/c 2 (area C) are fixed. 



In all fits, we fix the branching fraction of B° 



to the published value (33.2 ± 



3.0 ± 3.8) x 10 [2 If. Recently, Belle presented a preliminary result for the B 
branching fraction [22|, which is consistent with BaBar's value. We also float the B° 
a i =7r=F yield in the fit; the result with its error is consistent with the assumed value. 
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FIG. 1: Mi(7T7r)-M2(-7r7r) distributions in MC simulation of non-resonant 5° 
and B° — > p°p° decays, and the data, from left to right. Non-resonant B® 
MC decays are generated with 4- and 3-body phase space models. 



47T, B° -> p°7T7T, 

47T and B° — > p°iT7r 



The signal yields are extracted by performing extended unbinned maximum likelihood 
(ML) fits. In the fits, we use four dimensional (Mt> c , AE, Mi, M2) information to measure 
the branching fraction of B° decays into non-resonant 47r, p°p°, p°/o, /o/o and /ovr7r, and 
two dimensional (Mb c , AE) information for the decay B° — » p n7i. We perform a 2D 
fit to extract the p°7T7r yield since the Mi(7r7r)-M2(7T7r) region used in the fit is expected 
to consist of af 7r T , non-resonant An and p°7T7i decays and no contributions from other 
resonant decay modes. We define the likelihood function 

(\ Ncand / \ 
-e**) n E^ p i > w 
3 j »=i v 3 j 

where i is the event identifier, j indicates one of the event type categories for signals 
and backgrounds; rij denotes the yield of the j-th category, and Pj is the probabil- 
ity density function (PDF) for the j-th category. For the 4D fits, the PDFs are a 
product of two smoothed two-dimensional functions: Pj = Pj(M£ c , AE 1 , M{, MJ) = 
/smoothed (M^ c ,AE l ) x / smoot hed {M{ , M\) . For the B° -> p°mi branching fraction mea- 
surement, the PDFs are two-dimensional functions, i.e. Pj = Pj(M^ c ,AE l ) = 

/smoothed (M^ c , AE 1 ) . 

For the signal modes, the smoothed functions /smoothed (Mac; AE 1 ) and 
/smoothed (M[, M£) are obtained from MC simulations. For the M DC and AE PDFs, 
possible differences between the real data and the MC modeling are calibrated using a 
large control sample of B° — > D~ (Kh-k)^ decays. We find 16.2% (3.8%) of the recon- 
structed signal MC events are self-cross-feed (SCF) events for longitudinal (transverse) 
polarization; for SCF events at least one track from the decay B° — > p°p° is replaced by 
one from the accompanying B meson decay. We use different PDFs for the SCF events 
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and correctly reconstructed events, and employ the SCF fraction for a longitudinally 
polarized decay in the nominal fit. The uncertainty in the fraction of transversely 
polarized B° — > p°p° is taken into account as a systematic error. 

For the continuum and charm B decay backgrounds, we use a linear function for AE, 
an ARGUS function 23j for M^ c and a two-dimensional smoothed function for M1-M2. 
The parameters of the linear function and ARGUS function for the continuum events are 
floated in the fit. Other parameters and shape of the Mi-M 2 functions are obtained from 
MC simulations and fixed in the fit. 

For the charmless B decay backgrounds, we construct four separate PDFs for B° — > 
a^TT 11 , B + — > p + p°, B — > p + p~ + pn and other charmless B decays; all the PDFs are 
obtained using MC simulations. In the fit, while the yields of the B — ► afn T , B + -> p + p° 
and B — > p + p~ + pn are fixed to expected values obtained from measured branching 
fractions, the yield of other charmless B decays is floated. 



MEASUREMENTS OF BRANCHING FRACTIONS 

With the fitted signal yields n sig ., we calculate the branching fraction B using 

3=^-^ , (2) 

iV BB ' e MC " e PID 

where €mc is the overall reconstruction efficiency obtained using MC samples, epm is a 
PID efficiency correction that takes into account the efficiency difference between data 
and MC, and N B -g is the number of BB pairs. The production rates of B + B~ and B°B 
pairs are assumed to be equal. The PID efficiency correction is determined using an 
inclusive D* + — > D°tt + , D° — > K~tt + data sample where the track momenta and polar 
angles are required to be consistent with those of B° — > p°p°. 

The statistical significance is defined as S = a/— 2 ln(£ / £ max ) ? where Co and £ max 
are the likelihoods of the fits with the signal yield fixed at zero and at the fitted value, 
respectively. 

The 90% confidence level (C.L.) upper limit is calculated from the equation 



I N £{x)d 

XT c ( x ) d 



x 



X 



90%, (3) 



where x indicates likelihood variables corresponding to the yield, and iV is the upper 
bound for the yield that includes 90% of the integral of the likelihood function. 

The upper limit (UL) including systematic uncertainties is calculated by smearing the 
statistical likelihood function with a Gaussian, where the Gaussian width is the combina- 
tion of two total systematic errors: one is independent of the branching fraction and the 
other is proportional to it. The significance including systematic uncertainties is calcu- 
lated in the same way, but we only included the systematic errors related to signal yields 
in the convoluted Gaussian width. 

The fitted yields and branching fractions with systematic errors are listed in Table [B 
Fig. [2] shows the projections of the data onto AE, Mb c , Mi(tctt) and M 2 (7T7r) for non- 
resonant B° — > 47r decay in area A. Fig. [3] shows the projections of the data onto AE 
and Mb c for the B° — > p°nir decay in area B. Fig. H] shows the projections of the plots 
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for the B° — > p°p° decay in area C. We measure B° — > p°p° and non-resonant _B° — > 47r 
decays with 1.7a and 2.2a significance, respectively. There are no significant yields for 
B° - p°f , B° - /o/o, 5° - /ovttt or 5° - p^ix decays. 

TABLE I: Fit results for each decay mode listed in the first column. The signal yields, recon- 
struction efficiency, significance including systematic uncertainties (<S), branching fractions (B) 
and UL including systematic uncertainties are listed. For the yields and branching fractions, 
the first (second) error is statistical (systematic). 



Mode 


Yield 


Eff.(%) S B(xlO- 6 ) 


UL(xHT 6 ) (90% C.L.) 






Area A measurement: 




47T 


32.2 ± 14.91™ 


0.61 2.1 10.2 ±4.7+^ 


<17.3 






Area B measurement: 




p TT7T 


-11.5 ± 17.8^;° 


0.86 0.0 


<6.3 






Area C measurement: 




p°p° 

P°fo 


33.7 ± IQ.OtlH 
3.6 ±11.7±£f 


7.11 1.8 0.9±0.4t°;3 
4.01 0.3 0.2 ±0.6 ±0.4 


<1.6 
<1.0 


/o/o 


1.6 ±3.7 ±2.9 


2.23 0.4 0.1 ±0.3 ±0.3 


<0.8 


/ 7T7r 


-6.0 ± 19.7±H° 6 


0.71 0.0 


<8.6 



p° HELICITY 

We perform an angular analysis of B° — > p°p° using the sum of two p° helicity angle 
distributions for the signal candidates having M 1 (7T7r) and M 2 (7r7r) values in the signal 
region (0.626 GeV/c 2 < M h2 (nn) < 0.926 GeV/c 2 ). The p° helicity angle is defined as 
the angle between the 7r + direction and the B° direction in the p° rest frame. Fig. [5] shows 
the sum of the two p° helicity angle distributions; each bin has a B° yield obtained from 
a AE-Mb c fit that does not distinguish whether the B° decays into p°p°, afrc T or An. We 
fix the fractions of the three B° decays modes to values obtained from the the 4-D ML fit 
in area C. In the helicity angle fit, we vary the longitudinal polarization fraction for the 
B° — > p°p° component. A x 2 fit yields a longitudinal polarization fraction for B° — ► p°p° 
of Jl = 0.6 ± 0.2. The statistical error is obtained from a MC pseudo-experiment study, 
since the two p° helicity angle distributions are correlated. 

SYSTEMATIC ERROR 

The main systematic uncertainty for the branching fraction of non-resonant B° — > Att 
decay is the uncertainty of the B° — > afn T branching fraction. The main systematic 
uncertainties for the -B — > p°ittt branching fraction are the uncertainties of the non- 
resonant B° — > An and B° — > afn T branching fractions. For B° — > p°p° decay, the main 
sources of systematic uncertainties include the uncertainties on the branching fractions of 
B° -> af vr T , non-resonant B° -> 4vr, B° -> p°vr7r, 5 ± -> p°p ± and B° -> p°ir*°. 

We vary the branching fractions of 5° — > af 7r T (33.2 ± 4.8, in units of 10~ 6 ) [2l[ 
and B^ — > p°p ± (18.2 ± 3.0) [24J] by their ±la errors. On the other hand, we vary the 
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AE(GeV) M bl .(GeV/c 2 ) M,(jni) (GeV/c 2 ) M 2 (roi) (GeV/c 2 ) 



FIG. 2: Area A: projections of the four dimensional fit onto (a)AE, (b)Mb c , (c)Mi(7T7t) G (1.05, 
1.7) GeV/c 2 and (d)M 2 (vr7r) G (1.05, 1.7) GeV/c 2 . For the AE projection: 5.27 GeV/c 2 < 
M bc < 5.29 GeV/c 2 ; for the M hc projection: \AE\ < 0.05 GeV; for the M 1 ( 2 )(7T7r) projection: 
|A£| < 0.05 GeV and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 1.35 GeV/c 2 < M 2(1) (vr7r) < 
1.70 GeV/c 2 . The fit result is shown as the thick solid curve; the red solid curve represents the 
signal component, non-resonant B° — ► 47r decay; the blue dashed and and green dashed curves 
represent the cumulative background components from continuum plus b — > c backgrounds, and 
charmless B decays. 



(a) 



40 C 
30 
20 
10 

Rm -0.05 



Q 0.05 0.1 

AE (GeV) 




5.26 5.28 

M bc (GeV/c 2 ) 



FIG. 3: Area B: projections of the two dimensional fit onto (a)AE, (b)M bc . For the AE 
projection: 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 ; for the M bc projection: \AE\ < 0.05 GeV. 
The fit result is shown as the thick solid curve; the B° — > p°7r7r component cannot be seen 
here because of its negative yield; the blue dashed and and green dashed curves represent the 
cumulative background components from continuum plus b — > c backgrounds, and charmless B 
decays. 



branching fractions of B° — > p°K*°, non-resonant 5° — > 47r and 5° — > p°7T7r in the range 
(0,7.2) x 10~ 6 [24], (0,15.4) x 10~ 6 (area A measurement) and (0,4.1) x 10~ 6 (area B 
measurement), respectively, because these branching fractions are not measured with high 
significance. The fits are repeated and the differences between the results and the nominal 
fit values are taken as systematic errors. 

According to MC, the signal SCF fractions are 16.2% for B° -> p°p°, 10.6% for B° -> 
p°fo, 7.3% for B° -> f f , 12.4% for B° -> /ovrvr, 10.8% for 5° -> p°mi and 11.1% for 
non-resonant B° — > 47r. We estimate a systematic uncertainty for the signal SCF by 
setting its fraction to zero. A systematic error for the longitudinal polarization fraction 
of B° — > p°p° is obtained by changing the fraction from the nominal value fi = 1 to our 
measured value fi = 0.6. A MC study indicates that the fit biases are +7 events for 
B° -> p°p°, +1 event for B° -> p°/o, -1 event for 5° -> / /o, -3 events for B° -> /ovrvr 



10 




5.26 5.28 

M bc (GeV/c 2 ) 



M,(hti) (GeV/cT 



M 2 (i7i) (GeV/c') 



FIG. 4: Area C: projections of the four dimensional fit onto (a)AE 1 , (b)iWbc (c)Mi(7T7t) G (0.55, 
1.35) GeV/c 2 and (d)M 2 (7T7r) e (0.55, 1.35) GeV/c 2 . For the AE projection: 5.27 GeV/c 2 < 
M bc < 5.29 GeV/c 2 ; for the M hc projection: \AE\ < 0.05 GeV; for the M 1(2) (tttt) projection: 
\AE\ < 0.05 GeV and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 0.626 GeV/c 2 < M 2(1) (vr7r) < 
0.926 GeV/c 2 . The fit result is shown as the thick solid curve; the red solid line represents 
the signal component, B° — ► /5 /? decay; the blue dashed, dash-dotted and green dashed curves 
represent, respectively, the cumulative background components from continuum processes, b — > c 
decays, and charmless I? backgrounds. 




|cos9, 2 | 

FIG. 5: Sum of two p° helicity angle distributions for background-subtracted events in the region 
(0.626 GeV/c 2 < Mipfan) < 0.926 GeV/c 2 ). The points with error bars represent the B° yields 
obtained from AE-M^c ML fits for each helicity angle interval. The black overlaid curve is the 
sum of expected distributions from MC. The blue, cyan, red and magenta dashed curves represent 
the B° — > di7T, non-resonant B° — > 4-7T, longitudinally and transversely polarized B° —* p°p° 
decays, respectively. The red solid curve represents the sum of two polarized B° — > p°p° decays. 



and +6 events for non-resonant B — > An. We find that fit biases occur due to the 
correlations between the two sets of variables (AE, M\ )C ) and (Mi and M 2 ), which are not 
taken into account in our fit. We correct the yields in the fit for these biases and include 
the corrections as systematic errors. 

Systematic uncertainties for the AE-M^ PDFs used in the fit are estimated by per- 
forming the fits while varying the signal peak positions and resolutions by ±ler. System- 
atic uncertainties for the M\-M 2 PDFs are estimated in a similar way. We vary the mean 
and width of the p° and /o masses in the M\-M 2 PDFs for the decay modes B° — > p°p°, 
B° — > p°fo, B° — > fofo, and B° — > fann. For the non-resonant B° — > 4-7T, other rare B 



11 



decay, b — > c and continuum backgrounds, we vary their Mi-M 2 PDF shapes, repeat the 
fits, and take the deviations from the central values as the systematic errors. 



We test the possible interference between B° 



7r T , non-resonant B 



An and 



B — > p > p by toy MC. We add a simple interference model to the toy MC generation, 
which is, for p° — > 7r H 



m 2 — rriQ + imoT 



+ Ae 



Ti decay, modified from a relativistic Breit-Wigner function to 

2 

(m 2 — tuq) cos 5 — Tmo sin 5 



A Z + 2A 



(m 2 



m. 



I? 



+ {Tm f 



+ 



(m 2 



m. 



I? 



+ (iW 

(4) 

where A and 5 are the interference amplitude and phase, m and V are the p° mass 
and width, respectively. We assume that the interference term due to the amplitudes for 
B° — > af tt t and non-resonant B° — > An decays is constant in the B° — > p°p° signal region. 
Since the magnitude of the interfering amplitude and relative phase are not known, we 
uniformly vary these parameters and perform a fit in each case to measure the deviations 
from the incoherent case. The mean deviation is calculated, and we add and subtract 
the r.m.s. of the distribution of deviations from this value to obtain the systematic 
uncertainty. The interference systematic for B° — > p°p° decay is (^4 2 8 6 )%- 

The systematic errors for the efficiency arise from the tracking efficiency, particle identi- 
fication (PID) and TZ requirement. The systematic error due to the track finding efficiency 
is estimated to be 1.3% per track using partially reconstructed D* events. The systematic 
error due to the pion identification (PID) is 1.2% per track estimated using an inclusive 
D* control sample. The TZ requirement systematic error is determined from the efficiency 
difference between data and MC using a B° — > D + (K7i7i)ii~ control sample. Table II 
summarizes the sources of systematic uncertainties and their quadratic sum for each of 
the items. The overall relative systematic errors are (^39})% for the B° — > p°p° decay 
and (ii45)% for non-resonant B° — > An decay. 



SUMMARY 

In summary, we measure the branching fraction of B° — > p°p° to be (0.9 ± 0.4^q 4) x 
10~ 6 with 1.8a significance; the 90% confidence level upper limit including systematic 
uncertainties is B(B° — > p°p°) < 1.6 x 10~ 6 . Since no significant signal is found, we 
have assumed this mode is a longitudinally polarized decay (J'l = 1), to obtain the most 
conservative upper limit. For fi = 0.0, we obtain the central value of B(B° — > p°p°) = 
0.6 x 10~ 6 . Measurements of polarization and asymmetry in B° — > p°p° will be needed to 
improve the precision of the (f>2 constraint. 

On other hand, we find an excess in non- resonant B° — > Att decay with 2.1<r signif- 
icance. We measure the branching fraction and a 90% confidence level upper limit for 
non-resonant B° 4vr decay to be (10.2±4.7l 2 |) x 10" 6 and B(B° 4tt) < 17.3 x 10" 6 . 
This contribution was not taken into account in previous measurements of B° — > p°p° j^] . 
We find no significant signal for the decays B° — > p°fo, B° — > fofo, B° — > f ini and 
B° — * p°7T7r; the corresponding upper limits are listed in Table I. 
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TABLE II: Summary of systematic errors (%) for the branching fraction measurements, /l and 
fscF are the fractional uncertainties for longitudinal polarization and self-cross- feed. 



Source 


47T 


/9°7T7r 




P°fo 


/o/o 


/OTTTT 


Fitting PDF 


±6.8 


±9.6 


±13.1 ±197.7 ±159.4 ±194.4 


B(B° -► am) 
B(B° -► 4vr) 


±3.7 
— 


+100.0 
-96.5 
+181.7 
-86.1 


+4.7 
-5.3 
+14.5 
-7.1 


+2.8 
-5.6 
+16.7 
-52.8 


+18.8 
-12.5 
+37.5 
-50.0 


+48.3 
-51.7 
+326.7 
-146.7 


£(5° -► p 7T7r) 


— 


— 


-28.5 


-55.6 


±31.3 


-108.3 




— 


0.0 


±0.6 


±2.8 


0.0 


0.0 


fi(S u -> p [) K* [) ) 






±1.5 


±13.9 


0.0 


-16.7 


Sl 






-11.1 








fsCF 


-9.6 


±13.9 


-15.7 


-11.1 


±18.8 


±3.3 


Interference 






+12.6 
-4.8 








Fit bias correction 


+19.6 
-3.1 


±8.7 


+22.0 
-3.0 


+36.0 
-27.8 


±62.5 


+16.7 
-48.3 


Tracking 


±4.2 


±4.3 


±5.1 


±5.0 


±4.4 


±4.3 


PID 


±4.1 


±4.2 


±4.9 


±4.6 


±4.5 


±4.2 


1Z requirement 


±3.4 


±3.4 


±3.4 


±3.4 


±3.4 


±3.4 




±1.3 


±1.3 


±1.3 


±1.3 


±1.3 


±1.3 


Sum(%) 


+22.2 
-14.5 


+208.4 
-130.2 


+37.1 
-39.1 


+202.3 
-214.4 


+180.2 
-179.0 


+383.7 
-276.4 
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